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Abstract 

Using Langevin dynamics, we have investigated the dynamics of vortices in 
a disordered two dimensional superconductor subjected to a uniform driv- 
ing current. The results provide direct numerical evidence for a dynamical 
phase transition between a plastic flow regime and a moving "hexatic glass." 
The simulated current-voltage characteristics are in excellent agreement with 
recent transport measurements on amorphous Mo77Ge23 thin film supercon- 
ductors. 
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While the effects of disorder on the static Abrikosov lattice have been widely discussed 
|1]J2|], the driven disordered lattice has been much less studied, despite its obvious relevance 
to typical transport measurements. Schmid and Hauger || analyzed interactions between 
a pinning potential and the elastically deformed moving lattice. Several groups 0,[J have 
numerically demonstrated the importance of plastic deformation and associated dislocations 
in driven two-dimensional (2D) lattices. More recently, hysteretic, filamentary flow has been 
found numerically in the plastically deformed regime |J. 

An actual dynamic phase transition between the plastically deformed phase and a mov- 
ing lattice for the vortex lines was first proposed and numerically demonstrated in 2D by 
Koshelev and Vinokur JF|. Their idea is that the random "pinning noise" diminishes with 
increasing vortex velocity, allowing a high- velocity re-ordering. Because of the disorder, the 
moving lattice has been argued to be a novel "moving glass" phase with a finite transverse 
critical current M. Similar dynamical ordering phenomena may be relevant to Wigner elec- 



tron crystals ||, vortex lattices in Josephson arrays |nj, magnetic bubble arrays [|TTJ, and 
charge- density- wave systems JT2] , |TB| ]. 

Experimentally, a dynamical transition into a moving lattice phase has been suggested 
by transport measurements [Q, by neutron diffraction ]T5|] on flux lines in 2H-NbSe2, and 
by the current-voltage (IV) characteristics in the 2D amorphous superconductor Mo77Ge23 
||16|| . The latter displays maxima in dV/dl ||16|| , which may indicate a dynamical ordering 
transition. 

This Letter reports numerical studies of a model 2D vortex system which show evidence 
for three distinct dynamical phases. The evidence comes from calculations of both the IVs 
and the translational and hexatic order parameters. We also obtain a finite transverse critical 
current in the high-drive ordered phase, and a striking crossing of all the dV/dl curves for 
different temperatures at a single point (I cr , V cr ), in excellent agreement with experiment 

We consider N v two-dimensional vortices located at r-j in a thin superconducting layer 
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of thickness d, bulk penetration depth A(T) = A(0)/y 1 — T/T c , and fixed vortex density 
riB = I/as = B/(f)Q, where 0o — hc/2e. The classical Hamiltonian for the system is taken 
as 



^ = {E^(^f )+E^)}. (1) 



A±(T) 



The repulsive intervortex interaction is £/(x) = 8w 2x"t) 2 Kq{x), where Kq(x) is a summation 
of the modified Bessel function K (x) over image vortices under periodic boundary conditions 



in the xy-plane JT7j, and Aj_(T) = X (T)/d. The disorder potential V(r) is taken as N p 
potential wells of uniform depth U P (T,B) = a p (fi 2 ) /[64:7i 2 \±(T)] and radius r p . To carry out 
the calculation for a given pin configuration, we first anneal the system from T ~ T c down 
to T = lOOmK, in steps of AT = 200mK, using the Metropolis algorithm |Tj|. At each 
temperature T of interest, we take a snapshot J=0 of the vortex configuration during 
the annealing process, to be used as the initial dynamical configuration. 
The vortex dynamics are obtained from the overdamped equation 

rjr^t) = ff (t) + ff (f) + ^Jxi + f p (r,). (2) 

The first term on the right-hand side is the Brownian force due to Gaussian noise J7|. The 
remaining terms are the forces due respectively to the other vortices, the applied current, 
and an edge-smoothed version of the random pinning potential of eq. [I], as described in detail 
F^j . In actual runs, we used the Mo77Ge 2 3 [0 parameters A(0) = 7700 A, Ginzburg- 
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Landau parameter k = 140, d = 60 A (and hence A^ = 98.8/xm), T c = 5.63° K. Our bounding 
box containing 256 vortices for B = 0.5kG has an edge L x ~ 3.2/im, and we include up to 
200 x 230 image boxes within a circle, comprising a total linear dimension L ~ 6.5Aj_(0). We 
use a time step of 0.1 ~ lx t where t = — , an areal pin density B p j B = 4.0, and we 
take r p = 2£ ab (0), a p = 1. We also introduce a current density scale J p by y = U p (0)/r p . 
For our parameters, J p = 5.05 x 10 5 (A/cm 2 ), about 8 % of the depairing critical current 
density. Typically we discard the first 5000to of time steps to assure achievement of a steady 
state. 



To trace the disclination density rid(t) (i. e., sites with coordination number 7^ 6), we 
perform Delaunay triangulation on {ri(t)}^j every 10to- We also measure (typically aver- 
aged over r = 20000to) the vortex density-density correlation function G(r) = - J T dt < 
p(r,i)p(0,t) >, its Fourier transform 5(q), and the hexatic order parameter \I/ 6 [|19j. IV 
characteristics are obtained from the center of mass drift velocity of the entire vortex en- 
semble, and the Josephson relation E = n^0 o v / c x z. 

With no pins, we find that the lattice has a first order melting transition at T m = 0.41T C , 
with a Lindemann number of 0.19, and a simultaneous vanishing of both translational and 
hexatic order. This T m is reasonably close to the upper bound obtained from the dislocation 
unbinding picture: fcsT m = Ai(d0Q/128\/37r 3 A 2 (T m )), provided the renormalization constant 
A\ ~ 0.85 f2"0fl . With point pins, the vortices freeze into a disordered state with some residual 
hexatic order. An upper limit for the freezing temperature may be operationally defined 
as the point where the self-diffusion measure < | !■$(£) — r.;(0)| 2 > converges to a t-independent 
value. We find that the Tf ~ 2.0K ~ 0.87X m . 

Fig. [1] shows the simulated IV characteristics at several temperatures. The results at 
the lowest temperature (T/T m = 0.043) most clearly show that there are three distinct 
regimes J2TJ. (i) At low driving currents( J/ J p < 0.2), the lattice is pinned and the voltage 



is generated mainly by rearrangements of topological defects, which produce a very slow and 
unsteady advance of the pinned lattice. At least over the simulation time (~ 50, 000t ), the 
quenched-in hexatic order stays robust; we may therefore call this regime a pinned hexatic 
vortex glass (PHVG). (ii) At moderate J (0.2 < J/J p < 1.2), the IV's are strongly non- 
linear. The hexatic order is progressively destroyed with increasing J. The nonlinearity 
originates from the nature of plastic deformations in this regime, with highly filamentary 
vortex trajectories (cf. insets of Fig. |2|). Deep into this state, \l/ 6 ~ and the system is a 
plastic flow vortex liquid (PFVL). (iii) At high J, the lattice heals incompletely, and hexatic 
order is dramatically enhanced(\l/6 ~ 0.6) over the static phase (i) (^6 ~ 0.03). Translational 
order is still disrupted by a co-moving pattern of neutral disclination pairs. Both g(r) and 
S(k) display a strong anisotropy. We may thus call this state a moving hexatic vortex glass 



(MHVG). 

The simulated dV/dl curves shown in Figs. and |3] reproduce most of the essential 
experimental features |16|]. At the lowest T, Fig. |2| shows that the peak in the calculated 
dV/dl coincides with the dip in At higher temperatures the peak broadens and merges 
with the region at which \l/ 6 starts to increase again. In the simulations, the dynamic 
freezing transition (operationally defined as the sharp increase of ^§ at J ~ 0.8J P ) falls on 
the right-hand shoulder of the peak. In the experiment [fL6|l , the peak was interpreted to 



represent the increase in lattice correlation as it follows remarkably well the predicted form 
of Koshelev and Vinokur We therefore conclude that the peak in dV/dl marks the point 
where lattice correlations start to increase in both orientational and positional order. In the 
plastic regime [left shoulder of the peak], the density and tortuosity of the calculated vortex 
trajectories are very sensitive to both J and T (presumably via the T-dependent effective 
shear modulus of the vortex lattice). Qualitatively, the current paths resemble percolation 
paths just above the percolation threshold. We suggest that the current travels along an 
"infinite connected cluster" of open paths. These paths occupy an area which increases 
with increasing J, perhaps as a power law in (J- Jo), Jo being the current density where the 
infinite cluster first forms. 

As T — > T m , the dV/dl isotherms all tend to cross at a J cr ~ 0.12J P , strikingly reproduc- 
ing a prominent experimental feature ||16|| . The pinned fraction of the vortices also starts a 
sharp decrease at this same value of J / J p . Furthermore, as the inset of Fig. Q shows, this 
crossing coincides with the minimum in \l/ 6 . For J < J cr , we found that the vortex flow 
is composed of two components: unpinned individual vortices in the flux-flow state; and 
pinned vortex islands with local hexatic order. The latter contribute to the flow through 
intermittent motion. As J increases toward J cr , the size of the pinned islands shrinks and 
consequently, the moving vortex trajectories become less tortuous. We speculate that for 
each T, J cr (T) is the current density where the dV/dl for the pinned and unpinned con- 
tributions become comparable. The reasons for the crossing of the dV/dFs, and for the 
coincidence of this crossing with a minimum in ^q, will be discussed in greater detail else- 



where |22| . 



Finally, we consider the MHVG. In all our simulations, the moving healed lattice carries 
along bound disclination pairs ranging in density from 0.05 at T/T m = 0.043 to 0.4 at 
T/T m = 0.96. In the rest frame of the average velocity v, the lattice has local triangular 
order, but long-range translational order is disrupted by the disclination pairs (cf Fig. 4). 

This behavior suggests that the MHVG/PFVL is really dynamical freezing, as suggested 
by 0. To confirm this interpretation, we can estimate the effective "disorder temperature" 
Td- Viewed in the moving frame, the pins provide a random time-dependent potential 



which corresponds to an effective T^. For the present pinning potential, a simple analysis, 
including the vector nature of the forces, yields /c_eT rf = . Here a p = where 

n p is number density of pins. For our system, at zero temperature, U p = 7.5 x 10~ 15 erg, 
and a p = 1.02 x 10~ 5 cm. Using E = Bv/c, where E is the electric field, we finally obtain 
T d (K) « 6.6 x 10~ 12 /[^o(volt - sec/cm)]. For low temperatures (T/T m « 0.043), the 
PHVG/MHVG melting transition occurs near Et « 3 x 10~ 12 volt-sec/cm, or equivalently 
Td fa 2.2K, in excellent agreement with the simulated pin-free melting temperature. 

Note that in the MHVG phase, the lattice principal axis is not always parallel to the 
driven direction. The degree of misalignment in the annealed lattice correlates with the 
angle between the drift and driving directions || , suggesting the possibility of an anomalous 
Hall effect. This misalignment may occur because the lattice in MHVG is biased toward a 
direction which is frozen in during the J = annealing process, and thereafter frustrated 
from rotating parallel to the applied drive. To confirm this, we repeated with random initial 
configurations instead of an annealed state. The lattice then oriented its principal axis 
parallel to the driven direction, but the resulting state was higher in total energy than in 
the annealed case, with about twice the number of disclinations. 

Finally, we note the non-zero transverse critical current in MHVG. To see this effect, the 
lattice is first set in motion along a prescribed direction I. (of magnitude Je/J p = 2.6 in 
our simulations) for 20,OOOto- Then a small transverse current < Jt/Je < 0.01 is turned 
on along t = £ x z, and the lattice transverse velocity Vt is measured over 10, OOOto- The 

6 



results are shown in Fig. § for two angles 9 between I and a local principal axes of the moving 
lattice: 9 = and 9 = tt/6. For 6 = and T/T m = 0.043, v t drops abruptly at J t /J e = 0.002. 
For 9 = 7r/6, this drop is much less sharp. This transverse barrier seems to disappear very 
rapidly with increasing T, as shown for T/T m = 0.43, 9 = 0. A possible explanation for this 
behavior is suggested by the inset, which shows the quenched disclination patterns (open 
circle and squares) found in the moving phase. Presumably, during transverse lattice motion, 
these quenched defects slip, possibly reducing the transverse critical current below that 
expected from an elastic theory which does not allow such slippage. All these dynamic phases 
will probably manifest themselves more clearly in three dimensions. In a BiSrCaCuO-type 
materials, observation by neutron diffraction of a field-induced disordering transition at 
low temperatures was interpreted in terms of a topological phase transition from an ordered 



line lattice to a decoupled glass phase p4| . The glass phase has quasi-2D topological defects 
at low T, which may have IV similar to those described here. Thus, by driving the vortices 
with a high J at varying fields, and carrying out both resistivity and Hall measurements, 
one may be able to probe the moving lattice phase and its lower critical dimension. 

As this work was being completed, we became aware of work by Moon et al. on the 



similar model ||25|| . The authors gratefully acknowledge financial support by NSF Grant 
DMR94-02131 and the Midwest Superconductivity Consortium through DOE Grant DE- 
FG02-90ER-45427. 
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FIGURES 

FIG. 1. Simulated IV characteristics for T/T m = 0.043, 0.43, 0.65, 0.78, 0.87, 0.96, 1.087. 
Inset: J-T phase diagram constructed from the IV results and analyses of ^ (Fig. 2). The solid 
line (J/Jp = 1.2 + 1.2/(1 — T/T m )) and the broken lines are guides to the eye. Acronyms defined 
in the text. 

FIG. 2. Hexatic order parameter ^6 (broken line) and differential resistance dV/dl for 
T/T m = 0.043. Inset: vortex trajectories at points b and d in the plastic regime taken over 1000to- 
The dots represents those pinned vortices which remained static over this interval. Lower panel 
shows the time averaged g(r) function for J/J p : 0.04,0.4, and 1.98; they are marked a (PHVG), b 
(PFVL), and c (MHVG). 

FIG. 3. Simulated dV/dl for T/T m = 0.78,0.87,0.96, 1.087 and 1.30 near melting. Inset shows 
\&6- The crossing of all the dV/dl curves near J/J p ~ 0.12 coincides with the dip in \Ev 

FIG. 4. Transverse voltage induced by a transverse force z x J t || x. J t is turned on after a 
steady state is reached under the dominant force z x \\ y. Both T/T m = 0.043 (solid) and 0.43 
(broken) cases are shown. Inset: snapshot of vortices (points) for = 0,T/T m = 0.043. A pattern 
of disclinations with coordination numbers z = 8(thick circle), = 7(circle) and = 5(square) is also 
shown. 
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